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Abstract 

We present a model study of single-particle spectra and two-particle Bose- 
Einstein correlations for opaque sources. We study the transverse mass dependence 
of the correlation radii i?|| and i?o in the YKP parametrization and find a 
strong sensitivity of the temporal radius parameter to the source opacity. A 
simple comparison with the published data from 158 A GeV/c Pb+Pb collisions 
at CERN indicates that the pion source created in these collisions emits particles 
from the whole reaction volume and is not opaque. For opaque sources we find 
certain regions of inapplicability of the YKP parametrization which can be avoided 
by a slightly different parametrization for the correlator. The physical meaning of 
the modified parameters is briefly discussed. 



1 Introduction 



Bose-Einstein correlations in the two-particle momentum spectra of identical particle 
pairs provide a powerful tool to obtain information about the space-time structure of 
the particle emitting source. In a previous paper (referred to in the following as 
"Paper I") we began a detailed numerical model study of two-particle Bose-Einstein 
correlations for ultrarelativistic heavy ion collisions. We showed that collective expan- 
sion and temperature gradients lead to an M_|_ dependence of the correlation radii, 
and that by studying this dependence one can reconstruct the dynamical state of the 
collision fireball at freeze-out. 

A different feature of the particle source which can also cause observable effects 
in the M±_ dependence of the correlation radii but was not touched upon in Paper I 
is the source "opacity" ^ ^. This mechanism, which could affect the geometric and 
dynamical interpretation of the correlation radii, will be investigated here. 

"Opaque" sources emit particles from a thin shell near the fireball surface and are 
thus characterized by a smaller spatial extension of the emission zone in the "outward" 
(x) than in the "sideward" (y) direction^: (x^— y^) < 0. This is true in particular for 
emission functions from hydro dynamical simulations where freeze-out is implemented 
along a sharp hypersurface, characterized e.g. by a constant freeze-out temperature. 
Transparent sources of the type studied, for example, in P, ^ 0], on the other hand, 
feature a positive and generally small difference (x^ — y^) Q. 

In Refs. ||2|, ^ it was pointed out that this feature of opaque sources generally 
leads to smaller values for Rq than for Rg in the Cartesian parametrization of the 
correlator. In particular for pairs with vanishing transverse momentum K±^ = 0, one 
has i?^ = (y^) > R^ = (x^) if (5^ — y^) < 0. The existence of a thin emission 
layer with directed emission only into the outward hemisphere thus breaks the usual 
symmetry argument [P, |^ that \\tclk^^q{R^ — R^) = 0. That argument is based on 
the assumed azimuthal symmetry of the effective source for vanishing K± where the 
direction of the transverse pair momentum no longer serves to distinguish between 
the outward and sideward directions (parallel and perpendicular to K_|_); for opaque 
sources the orientation of the emitting surface itself provides that distinction. 

Here we show that in the Yano-Koonin-Podgoretskii (YKP) parametrization |lO| 
the opacity effects get enhanced in the "temporal" radius parameter Rq which turns 
negative for small K± and diverges to minus infinity in the limit K± ^ if (x^— y^) < 0. 
While this would destroy the interpretation of Rq in terms of an effective source lifetime, 
it would provide particularly clear evidence for surface dominated emission. Since no 
such evidence is seen in the data from the NA49 Collaboration [11, 12 1, rather stringent 
limits on the degree of "opaqueness" of the source created in these collisions can already 
now be established (see Sec. 3.3). 

In the course of this study we discovered that for opaque sources the YKP parame- 
trization may become ill-defined in certain kinematic regions. In Sec. |4.2| we introduce a 
modification of the YKP parametrization without this defect. The cost for the remedy 
is a less straightforward physical interpretation of its radius parameters. Fortunately, 
for transparent sources the YKP parametrization with its simpler space-time interpre- 
tation generally appears to work well. 

^The tilde notation is defined in Sec. ^ it indicates the variances of the distribution of emission 
points in space-time. 
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For the general formahsm and notation we refer the reader to Paper I. In the next 
Section we introduce the modification of the emission function from Paper I which is 
needed to parametrize the opacity of the source. In Sec. ^ we describe the results of 
our model study. Some problems with the YKP parametrization for opaque sources 
are discussed and corrected in Sec. and the results are summarized in Sec. |5[ As in 
Paper I we consider only directly emitted pions from the thermalized source, neglecting 
pions from resonance decays whose effects were studied in Refs. [|^]. 



2 A model for opaque sources 

Following the idea of Heiselberg and Vischer |||, ^ we introduce the opacity into the 
model emission function from Paper I via an additional exponential factor exp[— /ofj/A] 
which suppresses the emission of the particles from deep inside the source. lesii^,4') 
is the effective length which a particle emitted at point (r, ip) travels in outward (x) 
direction before leaving the source. We implement the Gaussian transverse density 
profile as follows: 
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A represents the specific mean free path of the particle in the medium. The actual mean 
free path would be obtained by dividing A by the medium density; this is, however, 
already included in Eq. (|l|). In the limit A — > oo the transparent source of Paper 
I is recovered. Note that, in contrast to Q, our mean free path A and transverse 
geometric radius R are time-independent. For later convenience we introduce the 
opacity parameter 
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Transparent sources are characterized by w = 0. 
The complete emission function now becomes 
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We do not specify the normalization of the emission function; in principle it is fixed 
by the total number of produced particles. For the shape of the single-particle spectra 
and two-particle correlations the normalization is irrelevant. As we will see a mean- 
ingful comparison with data is possible even without knowing the normalization of the 
emission function and permits us to exclude a large class of opaque source models. 

The collective flow velocity profile u{x) is parameterized as in Paper I. We will devi- 
ate from the discussion presented there by assuming a constant freeze-out temperature 
T. 
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3 One- and two-particle spectra from opaque sources 



In this Section we show and discuss the results of numerical calculations of correla- 
tion radii from the emission function (^). All calculations are based on the model- 
independent expressions which give the correlation radii in terms of space-time vari- 
ances of the emission function. For the Cartesian parametrization of the correlator 
they read 
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{Xfj_) and (• • •) denotes the space-time average taken with the emission 



function (see Paper I). For the YKP parametrization radii we have [0, |To| 
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In the following subsection we study the basic opacity effects on the one- and 



two-particle spectra. Some of these are analyzed in more detail in subsection 3.2. A 



comparison with available data is presented in subsection 3.3 



3.1 Basic opacity effects 

The basic free parameter of our study is the mean free path A. We will scan a rather 
wide range of opacities u) = R/X (u; = 0,1, 10). For r]f, which parametrizes the 
transverse flow at freeze-out linearly according to 

Vir) = Vf , (15) 

we will study two different values, r// = (a source without transverse flow) and 
r]f = 0.5 (rather strong transverse flow). All other model parameters are kept fixed in 
this and the following subsection at the values listed in Table [l| which are motivated 
by the SPS Pb+Pb data (see Sec. |3.3| ). The calculations in this and the following 
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Table 1: Model parameter values used in the calculations of Sec. 



3.1, 



temperature T 120 MeV 

average freeze-out proper time tq 8 fm/c 

mean proper emission duration At 2 fm/c 

geometric (Gaussian) transverse radius R 6.5 fm 

Gaussian width of the space-time rapidity profile Arj 1.3 

pion mass m^± 139 MeV/c"^ 



subsection are done for pions and pion pairs at midrapidity (Y^^ = Ytttt — Y^c 



0). 



To obtain an impression of how the source changes when varying lo and rjj we 
show in Fig. |l| for pions with transverse momentum K± = 400 MeV/c transverse 
cuts of the effective emission function, S{x,y;K± = 400MeV/c), as contour plots. 
When switching on the opaqueness, the "transverse" size (in y- or side-direction) of 



(a) 


_(b) 






(d) 


-(e) 




(f) 











20 
10 
I 
-10 
-20 
20 
10 
I 
-10 
-20 



-20 -10 10 20 -20 -10 10 20 -20 -10 
X (fm) X (fm) X (fm) 



10 20 



Figure 1: Transverse cuts of the emission function for midrapidity pions with transverse 
momentum Kj_ = 400 MeV/c (in x-direction) , for the model parameters given in 
Table |. Top row: no transverse flow (r]f = 0), bottom row: transverse How with 
r]f = 0.5. The columns correspond to different opacities: a; = (left), uj = 1 (middle), 
a; = 10 (right). 



the effective source is seen to increase dramatically (cf. Figs. |l|a-c); this is due to 
the suppression factor exp[— /gfr/A] which cuts out all of the interior of the source and 
leaves only the right hemisphere of the dilute tail of the Gaussian transverse density 
distribution. This cannot happen in the model of Heiselberg and Vischer |^ who use 
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a transverse box profile without Gaussian tails. For a source with the shape given 
in Fig. ||c the Gaussian approximation, on which the model-independent expressions 
from the beginning of this Section are based, may become questionable; for a qualitative 
understanding of the relevant features it should, however, be sufficient. 

Transverse flow (lower row in Fig. ||) is seen to decrease the effective source more 
in the sideward than in the outward direction; this effect is the stronger the larger the 
opacity w. 

The single-particle transverse mass spectra resulting from these models are shown 
in Fig. |2[ The normalization of the spectra is arbitrary since the emission function 




Figure 2: Comparison of spectra for transparent and opaque sources at midrapid- 
ity. For the model parameters see Table |l|. Thin lines correspond to models without 
transverse flow, thick lines to models with r]f = 0.5. For reference, solid lines repre- 
sent the spectra of the transparent model (to = 0). Dashed and dotted lines show the 
spectra from sources with uj = 1 and tj = 10 respectively. 



was not normalized. The interesting information is carried by the spectral slopes. For 
non-expanding sources, the slopes for opaque and transparent sources are identical. 
For expanding sources with fixed rjj, an opaque source yields much flatter spectra 
than a transparent source. The reason is obvious from Fig. ||: due to the opacity 
much stronger weight is given to the rapidly expanding surface than to the less rapidly 
expanding interior regions. This implies that the average flow velocity {v±) is larger 
for the opaque source. In it was argued that (for m± » rn-o) the slope of the 
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m_L-spectrum is given approximately by the blue-shifted temperature 



TeS 



1 - {vs_) 



(16) 



We should therefore expect roughly the same slopes if with increasing opacity lo the 



flow parameter rjf were reduced to keep {v±) fixed. In subsection 3.3 we will show how 
this works. 

The two-particle correlations resulting from these models are studied in Fig. ^. 
For midrapidity pions the cross-term Roi in the Cartesian parametrization and 
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Figure 3: Comparison of correlation radii from transparent and opaque sources for pion 
pairs with rapidity Y^^y^ = 0. The source parameters are given in Table ||. Thin lines: 
rjf = 0, thick lines: r]f = 0.5. Solid lines: uj = (transparent source), dashed lines: 



1, dotted lines: 



10. Note that Rs = R± and for y^^^ = the longitudinal 



correlation radius Ri = iiy and Y^^ = in LCMS. 



the Yano-Koonin velocity v of the YKP parametrization (^) vanish, and Ri = iiy. 
Furthermore, Rg = Rj_ in general. 

Surprisingly, without transverse flow the outward radius decreases only weakly with 
increasing opacity. This is a specific feature of our Gaussian parametrization of the 
source geometry and not true for a box profile as studied by Heiselberg and Vischer 
In their case no emission from regions outside the box radius R is possible, and 
opacity leads to an effective emission function which is "squeezed" into a very thin 
crescent-shaped region close to the edge of the box. In our parametrization, increasing 
opacity favors emission from the more dilute tail of the Gaussian density distribution, 
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and as a; increases the effective emission region just moves further outward, becoming 
part of a shell with a radius -Rsheii > R whose thickness happens to be more or less 
independent of to. The sideward homogeneity radius Rg, on the other hand, grows with 
i^sheii as seen in the top row of Fig. |l| and thus increases dramatically with increasing 
opacity. 

For transparent sources it is well known that transverse flow causes a reduction of 
the sideward and outward homogeneity regions (y^) and (x^), respectively, the more so 
the larger M± ^, 15, 17]. For fixed transverse flow parameter r]f, this effect is seen 



to be even stronger for opaque sources, reflecting again the growing average transverse 
flow velocity with increasing uj. The right upper panel of Fig. |3| shows that for very 
opaque sources this effect even eliminates the slight initial rise of Rq for small K±: the 
decrease of (x^) with transverse mass is then stronger than the increasing contribution 
from the term in @. 

It has been proposed in to identify opaque sources by looking for a negative 
difference R^ — R^, especially at K± ~ 0. In the lower left panel of Fig. ^ we show 
that much more dramatic opacity effects are seen in the "temporal" YKP parameter 
Rq. For midrapidity pion pairs Rq is given in the CMS frame (which is also the rest 
frame of the emitting fluid element where v(K) is zero) by 

Rl = {P)-^{xi) + ^{x'-y'). (17) 

For opaque sources the last term on the r.h.s. is negative and diverges to — oo as K± — > 
(/3j_ 0). Transverse flow reduces this negative contribution, but not sufficiently to 
reverse the overall negative sign of Rq in the region of small K± . As seen in Fig. ^ this 
specific feature of opaque sources is quite stable against a reduction of the opacity iv in 
the sense that it even shows for nearly transparent sources with a; = 1 at sufficiently 
large i^T^-values to be measurable. 

Compared to Rg, Ro, and Rq, the opacity effects on Ri = (shown in the lower 
right panel of Fig. ^) are weak. They can all be explained in terms of the different 
values for the average transverse flow velocities in the studied examples as discussed in 
Paper I. 



3.2 A more detailed discussion of — and R^ 

The interesting opacity effects on R^ — R^ and on Rq deserve a more detailed discussion. 
To this end we plot in Fig. Q the various contributions to these radius parameters from 
the space-time variances of a transparent (upper row) and an opaque source (lower 
row), both for the standard Cartesian (left column) and YKP parametrizations (right 
column). We show results for moderate transverse flow ly = 0.3 and midrapidity pion 
pairs, but we checked that no qualitative differences occur for forward rapidities. The 
diagrams in Fig. |^ can be concisely summarized by stating that for transparent sources 
the lifetime term (P) dominates the difference R^ — R^ and Rq (at least for small K±) 
while for opaque sources the difference (x^ — iP') of transverse spatial variances takes 
the leading role. 

These results confirm the conclusion of Ref. Q that, within the Cartesian parame- 
trization, opaque sources generically lead to negative values for the difference R^ — R^. 
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Figure 4: Opacity effects on — (left column) and Rq (right column). The solid 
lines show the radius parameters as they enter the correlation function, the other lines 
show how they are decomposed into contributions from various space-time variances. 
Left column: the dashed line shows (3']_{t^), the dotted line —2P±{xt), and the dash- 
dotted line (x^ —y"^)- In the right column the same symbols denote the same quantities 
divided by fi^. All curves are calculated at Y^^^ = 0, with model parameters from 
Table || except for i? = 7 fm and r/j = 0.3. Upper row: a; = (transparent source); 
lower row: a; = 10 (opaque source). 



Our calculations were done in the Yano-Koonin {v = 0) frameQ, and we checked that 
in this frame this qualitative conclusion is independent of the pair rapidity. In the 
YKP parametrization the same feature appears strongly enhanced in the parameter 
i?Q, especially at small values of K^, due to the division of the relevant combination 
of space-time variances by the factor j3\. A negative value for R^ corresponds to 
a negative value for R^ — R^ in the Yano-Koonin frame (which is usually not much 
different from the LCMS Q). 

3.3 Comparison with data 

In this subsection we compare the model calculations with recently published corre- 
lation data from 158 GeV/c Pb+Pb collisions at the SPS |11|, focusing on opacity 



^Remember that the Cartesian radius parameters Ro, Ri, Roi depend strongly on the longitudinal 
reference frame. 
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effects. To constrain the model parameters we make use also of preliminary data on 



single-particle spectra ||T2 



When we increase the opacity without changing any of the other model parameters 
we obtain a larger transverse HBT radius Rj_ since the opacity factor favors particle 
emission from the dilute tail of the transverse Gaussian density distribution. In order 
to reproduce a fixed measured value for R± with sources of increasing opacity we must 
therefore reduce the Gaussian width parameter R in the emission function. 

By the same mechanism larger opacities lead, at fixed rjf, to larger average trans- 
verse flow velocities in the source. These cause a flatter slope of the single-particle 
m_L-spectrum (Fig. ^) and a steeper M_L-dependence of the transverse HBT radius R± 
(Fig. ^). On a qualitative level both effects can be studied analytically: for pions the 
inverse slope parameter of the single-particle m_L-spectrum is given by Eq. (p^), while 
the M_L-dependence of R± is approximately given by 



15[ 
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R^ 



1 + 



(^) K)^ 



(18) 



This expression is derived from Eq. ( pT|) by evaluating the average (y^) over the emission 
function via saddle point integration 15 1. Although not quantitatively reliable 
it illustrates nicely an important point 
each by itself, ambiguous and do not a 



1|, Ig, P, |n|, |12|: while (H) and (|18|) are, 
low to separate T and {v±), the correlation 
between those two parameters is opposite in the two equations. By combining them the 
collective transverse flow {v±_) can be isolated from the random thermal motion ~ T. 
The value for T given in Table |l| stems from a roughQ numerical fit |21] to the measured 
single-particle spectra and correlation radii from Pb-|-Pb collisions at the SPS 18|. 
We expect it to be accurate to about ±10 MeV. 

Once T is known, the average transverse flow is fixed by the slope of the single-pion 
transverse mass spectrum. For T = 120 MeV the NA49 spectra |l^ yield an average 
transverse flow velocity {v±) ~ 0.4c — 0.45c [p2| , pO| , [2l[| . For our source model the 
average transverse flow velocity is given by (^ = r/R) 
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This equation implies that, for fixed {v±), rfj must be reduced if oj is increased. There- 
fore, if we want to keep the absolute value of R±, its slope, and the slope of the 
single-particle m_L-spectrum fixed, we must decrease both R and r]f as we increase the 
opacity oj. 



^In the fit resonance decay contributions were only included in the single-particle spectra, but 
not in the calculation of the correlation radii. Contrary to the analysis in |^[ we evaluated, however, 
the spectra and HBT radii numerically with the source (^) (for lu — Q) instead of using the simple 
approximations ( |l^ and (p^). 
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Table 2: Transverse flow strength r/j and Gaussian radius R for the three models which 
are used in the comparison with the data in Sec. 3.3. 



uj (opacity) 




R 


(transparent) 


0.4 


5.74 


fm 


1 


0.345 


4.83 


fm 


10 (opaque) 


0.215 


3.35 


fm 



Three models with opacities u; = 0, 1, 10 were studied and compared with the data. 
In each case we tuned r]f such that {v±) = 0.44 was kept fixed, and we adjusted R 
such that the size of R±{M±) is roughly reproduced. This choice was motivated by the 
NA49 data shown in the right column of Fig. 4 of Ref. [0] (HBT radii from 158 GeV/c 
Pb+Pb collisions, h^h^ correlations in the rapidity window 3.9 < 1^,^ < 4.4). The 
values of R and ijf used in these calculations are listed in Table |2[ To fit the parallel 
correlation radius (see later) the average freeze-out time tq has been set to 6.1 fm. 
Other model parameters are taken from Table [l|. 

In Fig. |5|we plot the single-particle m_L-spectra resulting from these models. (Recall 




Figure 5: Single-particle m_|_-spectra from the models listed in Table ^. Solid, 
dashed, and dotted lines correspond to models with uj = 0, 1, 10, respectively. 
The thin dash-dotted auxiliary lines show a straight exponential with inverse slope 
Teff = T^{1 + (?;_l))/(1 - {v_l)), for T = 120 MeV and {v_l) = 0.44 resp. {v_l) = 0.49. 
The spectra are computed at midrapidity where the analysis of was performed. 
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that only the slopes, but not the normalization are relevant here.) According to Eq. (16) 
we expect the same slopes for all spectra, which should be characterized by the blue- 
shifted temperature with {v±) = 0.44. This is roughly borne out by the calculation, 
but for transparent and mildly opaque sources the spectra are visibly concave |16| and 
show an increasingly flatter slope for m± > 1 GeV. The reason for this is well-known 



16]: for large m± the spectra are dominated not by particles from fluid cells with the 



average transverse flow velocity {v±), but from more rapidly moving fluid cells in the 
tail of the transverse density profile. This detail cannot be accurately accounted for by 
the simple formula (|l^). In other words: for a specified p± there is only a part of the 
source contributing to the particle production - this is the effective source for that p±. 
The slope of the spectra at that value of p± is given by the blue-shifted temperature 
where the average transverse expansion velocity is to be calculated over the effective 
source. For the very opaque source this average velocity is (accidentally) independent 
of p_L, and the spectrum is thus exactly exponential. However, to reproduce its slope 
via ( p^ ) we need {v±) = 0.49, a somewhat larger value than the average transverse 
expansion velocity computed over the complete source. Note that the experimental 
spectra resemble more the concave spectrum of the transparent source than the straight 
line originating from the opaque one. 

Fig. ^ shows the corresponding set of YKP correlation radii, evaluated at = 
1.25 where the NA49 collaboration has published the M± dependence of the three 
YKP radius parameters from 158 A GeV/c Pb+Pb collisions at the SPS (Fig. 4 in 
Ref. |l^). Note that the systematic error of the data, estimated to be ±15% |11|, 
is not included in the plots. In Fig. ^ one sees that by simply adjusting the source 
parameters to keep the magnitudes of R± and {v±) approximately constant, the M±- 
dependence of R± is reasonably well reproduced by all three models. This implies that 
the average transverse flow velocity is indeed the dominating factor for the M_L-slope 
of the transverse HBT radius. We found that the agreement between model and data 
can be further improved by additional fine-tuning of the parameters, in particular by 
taking into account a non-vanishing transverse temperature gradient. Since a fully 
quantitative analysis requires, however, also the inclusion of resonance decays which 



may cause an additional M± dependence [13, 23 1, and since the following arguments 



will not depend on such details we will not enter into this discussion here. 

The longitudinal radius parameter shows even less sensitivity to the opacity 
uj (Fig. |6|b). The reason is that it is determined mainly by the values tq (average 
freeze-out proper time) and At (mean emission duration), and its M± dependence is 
dominated by the boost-invariant longitudinal expansion p^ , |^. Since, except for very 
forward and backward rapidities, opacity affects the emission function mostly in the 
transverse direction, is hardly modified. By appropriately adjusting tq, the data 
for are very well reproduced. Note that for opaque sources the calculated curves 
feature a gap at small values of M±_ and behave strangely near M± = rriT^; the origin of 
this behaviour will be explained in Sec. ^. This problem does not arise at midrapidity 
where unpublished data on the M_|_-dependence of can be found in the PhD theses 
1 25] which we were able to fit equally well all the way down to K± = with our emission 
function for all three values of uj, using tq = 8 fm/c. 

Since the space-time rapidity width Arj = 1.3 is fixed by the width of the single- 
particle rapidity distribution ]|l^, |ll| and tq has been set to the value of 6.3 fm, at this 
point the emission duration At is the only remaining free parameter in the emission 
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Figure 6: M±_ dependences of the YKP correlation radii at Y^y^ = 1.25 resulting from 
the three models in Table ^. (Same line symbols as in Fig. |.) The dat a are from 158 
A GeV/c Pb+Pb collisions taken by the NA49 collaboration pl[|. 



function. For transparent sources (w = 0) it should be fixed by the asymptotic value of 
Rq at large values of M_l ||1^. (Its effects on i?_L are negligible, and for Ar <C tq they 
are also small for Q.) In Fig. |6|c we compare with the data our results for i?Q for a 
fixed value Ar = 2 fm/c, but for three different opacities. Clear disagreement with the 
data is seen for the opaque models. Already the still rather transparent model with 
uj = 1 (cf. Fig. HI) misses nearly all data points. 

The question arises to which extent this disagreement can be avoided by readjusting 
the duration of particle emission Ar. As mentioned before and seen in Fig. for 
opaque sources the asymptotic value of Rq at large M_\_ is no longer given by the 
effective source lifetime (P) in the YK frame alone, but also receives a strong negative 
geometric contribution {Sp' — y^) (cf. Eq. (IT^)). One might try to compensate this effect 
by increasing for opaque sources the lifetime Ar. This is physically not unreasonable 
because opaque sources emit particles only from the surface; to produce the same 
total particle yield the resulting reduced brightness must be compensated by a larger 
emission duration. However, larger values of Ar lead also to an increase of R\\ which 
must be compensated by reducing the average freeze-out time tq. Since our model 
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parametrization looses its meaning for At > tq, there are clear limits to what can be 
achieved in this way. 

In Fig. 1^ we show and Rq for a mildly opaque source with u = 1 for different 




Figure 7: The YKP radius parameters Rq (a) and (b) computed for a; = 1 and 
different values of tq and At at Y^^^ = 1.25. Solid lines: tq = 5fm/c, At = 3fm/c; 
dashed lines: tq = 5fm/c, At = 4fm/c; dotted lines: tq = 6fm/c, At = 3fm/c; dash- 
dotted lines: tq = 7fm/c, At = 3fm/c. Other model parameters are taken from the 
corresponding calculation in Fig. ^. 



combinations of tq and At and compare them with the data [11|. For R and r/j we 
took the corresponding values from Table |2[ 

One sees that the largest lifetime At compatible with the laTge-K± data for is 
4 fm/c. To maintain a reasonable description of in this case tq must be reduced 
to 5 fm/c which brings it dangerously close to At. Still, the corresponding curves for 
Rq all miss the first data points by a large margin. Since however, as already seen in 
Fig. ^, at the given value of the pair rapidity the calculated curves for Rq feature a gap 
at low K± (indicating that for our model the YKP parametrization is not well-defined 
in this momentum region, see Sec. one may question the significance of this failure 
of our model to reproduce the data. For this reason we also checked the behaviour of 
Rq and in different rapidity bins, using the unpublished data of Ref. [25|. In all 
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cases the data points for Rq at the lowest K± values are positive, while the calculated 
curves for opaque sources with uj > 1 give strongly negative values, missing the data by 
far. In particular, this is true at midrapidity where the YKP parametrization remains 
well-defined down to K± = even for strongly opaque sources. 

We conclude that the NA49 data appear to exclude source opacities a; > 1, i.e. 
sources with a "surface thickness" X < R. The pion source created in Pb+Pb collisions 
at the CERN SPS seems to be rather "transparent" at freeze-out, i.e. at the end of the 
hydrodynamical expansion phase the pions freeze-out in bulk, decoupling essentially 
everywhere at the same time, including the center of the collision region. It must be 
kept in mind, however, that this conclusion rests heavily on the positive experimental 
values for R^ at small K± reported by the NA49 collaboration. In ||2|, |3| it was remarked 
that preliminary NA44 data from Pb-|-Pb collisions at the SPS indicate that at small 
M_i_ in the Cartesian parametrization R^ might be smaller than R^ (see |26]); this would 
favor an opaque source. In our opinion it is premature to draw firm conclusions from 
these data because, unlike NA49 ||ll[, the NA44 Collaboration has not yet presented 
a YKP fit of their data, where the opacity signal should be much clearer. Taking 
their Cartesian fit presented in [26| at face value we conclude that, while a mildly 
opaque source with a; ~ 1 cannot be excluded, the NA44 data are also inconsistent 
with strongly surface dominated emission {uj > 10). 

In the next Section we will explain the origin of the missing pieces in the curves for 
a; > 1 in Figs. |6| and 0. 



4 Modified YKP parametrization 

In the previous Section we discovered an intrinsic deficiency of the YKP parametriza- 
tion which had gone unnoticed before. In this Section we explain the limitations of the 
YKP parametrization and suggest a slightly different parametrization which avoids 
these problems and should be particularly suitable for strongly opaque, azimuthally 
symmetric sources with dominant longitudinal expansion. 



4.1 Forward rapidity pion pairs from opaque sources 

In Fig. ^ we show the YKP radius parameters Rq and R\\ and the YK rapidity Y^^ = 
^ ln[(l-|-T;)/(l— w)] (relative to the LCMS) for pion pairs at forward rapidity Y^^^ = 1.5 
for opaque sources with and without transverse flow. Let us first discuss the Yano- 
Koonin rapidity Y-y-^^ : For transparent sources we know from [0, ^] that the YK rapidity 
is negative in the LCMS, i.e. that the effective source moves somewhat more slowly in 
the beam direction than the emitted pairs. For opaque sources this seems no longer 
to be true; for example, in Fig. |^ in the absence of transverse flow the YK rapidity 
comes out positive in the LCMS. We hasten to stress that this does not imply that in 
this case the effective source moves faster than the emitted pairs (which would indeed 
be counterintuitive); it rather reflects the fact that for opaque sources the geometric 
correction terms in the expression for the YK velocity (see Eq. (4.1) in Ref. Q), in 
particular the one proportional to (x^ — y^), are large and spoil the interpretation of v 
as the longitudinal source velocity. 

Figs. ^ and |8| show that these correction terms can cause even more severe problems: 
in certain Mj^-regions the argument of the square root in (^) turns negative and the YK 
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Figure 8: The YKP radii and the Yano-Koonin rapidity in the LCMS for pion pairs 
with y^j^ = 1.5 from an opaque source (a; = 10) with = 4 fm. Sohd hnes: r]f = 0. 
Dashed hnes: r/f = 0.6. The other model parameters are taken from Table ||. The 
arrows connect the curves across the gap in which the YKP parametrization is not 
defined. 



velocity (and thus also Rq and see (^), (|To|)) becomes undefined. The conditions 
for the occurrence of such a pathological behaviour are discussed in Appendix It 
appears that it is connected in a generic way to (x^) < (y^) which, as we have seen, 
characterizes opaque sources. 

Note that the mentioned "strange" behaviour of the YK velocity, in the regions 
where it is defined, is reflected via Eqs. @ and (|l0|) in also unexpected behaviour of 

and Rq especially near to the gap. 

4.2 A solution of the problem 

This problem of the YKP parametrization can be avoided by a simple technical modifi- 
cation which we call "Modified Yano-Koonin-Podgoretskii" parametrization. It differs 
from the YKP parametrization by using rather than qj_ = \/ + qf, as one of the 
three independent components of the relative momentum q. For distinction we denote 
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the corresponding HBT parameters with a prime: 

C7(q, K) = 1 + exp[-ii^ (K) ql - (K) (g? _ (gO)2 
-(<(K)+i?f(K)) {q-U'm 

with 

C/'(K)=7'(K)(1,0,0,^;'(K)). 
Introducing in analogy to ([l2|)-(|l4|) the shorthands 

2 1 



A' 
B' 
C 



P. 



-{tx) + 

Pi 



f3l 



{z')-2^{zx) + ^{x'), 



PL 
PI 



1 



{zt)-^{tx)-—{z±) + ^{x') 

PI. P± P± 



(20) 

(21) 

(22) 
(23) 
(24) 



the modified YKP parameters can be expressed by the same formulae as the original 
ones: 





A' + B' 1 


v' 




2C' y 




= B'-v'C 


Uq 


= A!-v'C' 


R'l 


= 



A' + B' 



(25) 

(26) 

(27) 
(28) 

In the modified YK frame, defined by v' = 0, we thus have R'^ = B' and Rq = A'. 
Of course, since both YKP and Modified YKP parametrizations are just two ways of 
parameterizing the same correlation function, they are related in a simple way: 



A 
B 
C 



1 Rq +1 V R\\ - -^R 

P± 



1 



J2 



12 T>I2 , J2j2 p/2 Pf p/2 

P_L 

[R'l + R'^)^%'-§R'l. 

Pi 



The inverse relations are given by 

A' 

B' 

C 



-i'^Rl + Yv'Rt + ^R 



2„,2 d2 



1 



-f^Rl + -f^^Rl + 



PI 

''rI 

Pi u2 



PL 

PI 



{Rl+Rl)^\ + ^Ri. 

Pl_ 



(29) 
(30) 
(31) 

(32) 
(33) 
(34) 



In Appendix ^ we show that the Modified YKP parametrization is defined ev- 
erywhere except for the point K± = to which it can be smoothly extrapolated. 
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Furthermore, whereas the relative momentum components used in the original YKP 
parametrization satisfy the inequality 

q±>qo = -3-9 - -^Qi , (35) 

which means that the data points never fill the whole three-dimensional g-space, no 
such restriction exists for Qs which is used in the modified parametrization. This should 
help to avoid certain technical problems in the fitting procedure which can occur with 
the YKP parametrization j28| j. In any case, it may be useful to check the YKP fit 
against a modified YKP fit via the relations (p9|)-(|3T|), in order to avoid the pitfalls 
related to the possible non-existence of a YKP parametrization for the data under study 
(which may not show up clearly in the fitting process but might cause it to converge 
to a wrong result). We would strongly recommend this check in order to support 
our conclusion from the previous section that the data exclude opaque sources; this 
conclusion was based on the value of Rq in a region where, if the source were indeed 
opaque, the existence of YKP parametrization might be questionable. 

Unfortunately, the physical interpretation of the modified YKP radii is no more 
as straightforward as that of the original ones. Wherever in A, B, C the difference 
(x^ — y^) occurs it is now replaced in A' , B' , C by (5^) alone. For transparent sources 
the difference (x^ — y"^) is usually small and the corresponding correction terms 
to the leading geometric contributions to Bjq and BI are of minor importance; this 
permits a direct interpretation of Rq and R\\ as the effective lifetime and longitudinal 
size of the source in its own rest frame For such sources the occurrence of {Sp') 
alone in the modified radius parameters is certainly a drawback and usually leads to 
large corrections which invalidate a naive geometrical space-time interpretation (see 
Fig. ^). For opaque sources the appearance of (x^) (which is related to the curvature 
and thickness of the emitting surface shell) instead of the generically much larger 
combination (x^ — y^) may at first sight appear as an advantage. Nevertheless, unless 
the emitting surface layer is indeed very thin (and very flat!), the resulting correction 
term in particular to the leading contribution in R'q cannot be neglected and spoils its 
simple interpretation as an effective source lifetime. This is shown in the left upper 
panel of Fig. P; even for opaque sources with u) = 10, one sees that R'q is in most of the 
cases much bigger than the effective source lifetime Ar = 2 fm/c, even at large values 
of Mj_ (where /3± 1). 

The effects on the modified longitudinal radius parameter (right upper panel in 
Fig. ^ are much smaller, due to the small value of the longitudinal pair velocity f3i in 
the modified YK frame (where v' = 0) which multiplies the correction terms. (5i is small 
in the modified YK frame because, like the original YK rapidity, rises linearly with 
the pair rapidity Y^j^ with nearly unit slope, reflecting the boost-invariant longitudinal 



expansion of the source (see upper row of Fig. 10). The difference between Y^^^ and 
the longitudinal flow rapidity l^gpg around the point of maximum emissivity (see lower 



row of Fig. 10) is somewhat larger than for the original YK rapidity for transparent 
sources [|10|, |7[, but still small enough to consider as a good approximation for the 
rapidity of the effective source. The difference between Y^k ^^'^ rapidities of the 
LCMS and LSPS disappears in the limit K±^ — > cxd, as for the original YK rapidity. 
Opacity effects on Y^^ are seen to be small. 
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Mj^ (MeV/c) Mj^ (MeV/c) 

Figure 9: The correlation radii R'q^ and of the Modified YKP parametrization, 
for pions at Y^^^ =1.5. Dashed and dash-dotted hnes correspond to the transparent 
source without transverse flow (dashed) and with rjj = 0.6 (dash-dotted). SoHd and 
dotted lines correspond to an opaque source (w = 10) with rjf = (sohd) and r]f = 0.6 
(dotted). The lower row shows the modified Yano-Koonin rapidity Y^^ in the LCMS 
(left) and relative to the LSPS (right) as a function of M±. Other source parameters 
as in Fig. 



5 Conclusions 

The most important result found here is that a significant opaqueness of the source leads 
to dramatic effects on the YKP radius parameter Rq. Opaqueness was parametrized 
by the opacity, i.e. the ratio between the Gaussian transverse geometric radius R and 
the surface thickness or mean free path, uj = R/X. Even for moderate values of uj, 
u; ~ 1, i?o becomes negative at small, but experimentally easily accessible values of 
Kj_. For higher opacities this feature extends over larger M± regions. For u; > 1 the 
effect should be clearly visible in the existing data from Pb-|-Pb collisions at CERN. 

However, we also encountered a problem which made the comparison with data 
slightly problematic: we discovered previously unknown regions of inapplicability of 
the YKP parametrization of the correlation function. These regions of inapplicability 
are connected with the use of q± = y/q'i + q'o instead of Qs as an independent relative 
momentum variable. We found that for opaque sources the YK velocity may become 
ill-defined. For this situation we suggested a modified YKP parametrization which is 
always well-defined and particularly suited for opaque sources. Even for transparent 
sources it can always be used as a technical tool to check the correct convergence of 
the YKP fit, by using the relations (p9[)-(^). On the physical level, the modified YK 
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Figure 10: The modified YK rapidity Y^^ in the CMS (upper row) as a function of the 
pion pair rapidity in the CMS, for two values of K± = 150 MeV/ c (left) and 550 MeV/c 
(right). Line symbols and model parameters as in Fig. 



velocity v' can again be interpreted, in good approximation, as the longitudinal velocity 
of the effective source. It continues to reflect the longitudinal expansion of the source 
through a strong dependence on the rapidity of the emitted pairs. The interpretation 
of the modified transverse and longitudinal radius parameters R'j_ and i^y as transverse 
and longitudinal regions of homogeneity in the source rest frame remains valid with 
sufficient approximation. The interpretation of Rq (which is now always positive) as 
the square of the effective source lifetime, however, is spoiled by a large geometric 
correction. The extraction of a reliable estimate of the source lifetime for opaque 
sources thus appears to be a very difficult problem. 

From the fact that the published data |jll|] show only positive values for Rq we 
concluded that the experiment favors a source with volume-dominated emission. Pion 
freeze-out in heavy-ion collisions thus appears to be similar to the decoupling of the 
microwave background radiation in the Early Universe: at a certain point in time, 
when the matter has become sufficiently dilute and cool, suddenly the entire fireball 
(universe) becomes transparent. However, since for large opacities and in the rapid- 
ity region covered by the published data |11] the YKP parametrization was found to 
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break down in the critical region of negative Rq, two additional checks are required to 
confirm this conclusion: the M^-dependence of the YKP radii at midrapidity should 
be published, and a careful cross-check of the fitted correlation radii with the Carte- 
sian and/or modified YKP parametrizations should be performed, using the relations 
given in Sec. |^ and in ||^, |l^. In addition to boosting our confidence in the correct 
convergence of the multidimensional fit to the measured correlation function, such a 
cross-check would also exclude possible doubts about the applicability of the YKP 
parametrization to the data. Without doing that one can still try to draw conclusions 
based on the difference — R^ at small K± , but compared to the opacity signal carried 
by Rq this difference is expected to be smaller and more difficult to measure accurately. 
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A Definition range for the YKP parameters 

The Yano-Koonin velocity v in (P) is only defined if the discriminant 



D = {A + Bf- 4^2 



(36) 



is positive. Here we study the conditions under which this is the case. To this end we 
write the expressions (12)- (14) as 



A 
B 
C 



A' -'(f), 



B' 



If 
PI 



' if) 



(37) 
(38) 
(39) 



We then have 



and thus 



We will prove that 



\A' + B'\ =A' + B'. 
{A' + B' f - 4C"^ > 



(40) 

(41) 
(42) 

(43) 
(44) 
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or 



\A' + B'\ -2\C'\ > 0. 
Due to (1^) the inequality (|45|) can be written as 

A' + B' ^ 2C' > 



(45) 



(46) 



where the upper sign stands for C" > and the lower for C' < 0. Inserting expressions 
(|40|)-(E^) we obtain 



(47) 



This proves (0). Since A', B', and C are the shorthands belonging to the modi- 
fied YKP parametrization, we have also proven that this parametrization is defined 
everywhere (except, of course, the point = /?± = 0, 0]). 
Prom expressions (B^-(B9|) we see that 



A + B = A' + B 



2C = 2C'-^{y' 



PI 



It is therefore possible to obtain negative values for D: 

D = {A + Bf - AC'^ <{). 



(48) 
(49) 

(50) 



For example, if /3/ = and {y ) is fixed at a value bigger than (x ), one can find values 
for (3^ such that 



\A + B\ 



< 2\C'\ = 2\C\ 



(51) 



We must therefore conclude that there are kinematical regions where the YKP para- 
metrization is not defined. 
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